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The mass-dependent symmetry energy coefficients a„ m (A) has been extracted by analysing the heavy nuclear mass differences 
reducing the uncertainties as far as possible in our previous work. Taking advantage of the obtained symmetry energy coefficient 
a„m(A ) and the density profiles obtained by switching off the Coulomb interaction in 208 Pb, we calculated the slope parameter L o n 
of the symmetry energy at the density of 0.11 fm '. The calculated Lo.n ranges from 40.5 MeV to 60.3 MeV. The slope parameter 
Lo.ii of the symmetry energy at the density of 0.1 I fm is also calculated directly with Skyrme interactions for nuclear matter and 
is found to have a fine linear relation with the neutron skin thickness of 208 Pb, which is the difference of the neutron and proton rms 
radii of the nucleus. With the linear relation the neutron skin thickness A R np of 208 Pb is predicted to be 0.15 - 0.21 fm. 
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1 Introduction 

Symmetry energy is one of the basic features of the equation 
of state (EOS) of nuclear matter. It represents the energy 
cost in translating all the protons to neutrons in symmetric 
nuclear matter for per nucleon approximately. The density- 
dependent symmetry energy S ( p ) has attracted great attention 
in nuclear physics and astrophysics such as heavy ion reac¬ 
tion EMI, the stability of superheavy nuclei im nuclear 
structure mm, the structures, composition and cooling of 
neutron stars 1T&142T11 and even some new physics beyond 
the standard model |22][23]. Consequently, the slope and 
curvature parameters which describe the density-dependence 
are of much importance to understand a variety of issues in 
the mentioned areas. The slope parameter L is vital partic¬ 
ularly, which mostly decides its density-dependence. At the 
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saturation density, the slope parameter L probed by plenty of 
methods can vary largely 1241 . Fortunately, recent available 
constraints on L from terrestrial laboratory measurements and 
astrophysical observations are in agreement with L - 55 ± 25 
MeV, which is summarized by Chen recently l25l . But 
further efforts are still needed in the determination of the 
density-dependence of the symmetry energy, in particular, 
the high density behaviors. Although the nuclear properties 
given by various effective interactions are very different, the 
finite nuclei structure characters provided by them are ap¬ 
proximately unanimous. For this reason, it is a significant 
way that we explore the nuclear matter with the help of prop¬ 
erties of finite nuclei. In our work, we utilize the properties 
of finite nuclei to probe the density-dependence of symmetry 
energy around the subsaturation density of 0.1 lfnT 3 . 

In recent years, it has been established that the neutron skin 
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thickness A R np of 208 Pb is linearly correlated with the density 
dependence of the nuclear symmetry energy around satura¬ 
tion 1261 - 1311 . In the present work, we also observe that there 
is a strong linear relation between the neutron skin thickness 
A R np of 208 Pb and the slope parameter Lo.n of the symme¬ 
try energy at the density of 0.1 lfnT 3 . According to this we 
can obtain a strong constraint on the density dependence 
of S (p) with a measurement of A R np with a high accuracy. 
Correspondingly, the neutron skin A R np can be predicted by 
the the slope parameter L^’j l . As the smaller neutron skins 
in heavy nuclei tend to yield smaller neutron star radii {32) 
given by Horowitz and Piekarewicz and the value of neutron 
skin play a key role in deciding whether the 1,4M Q neutron 
stars can have a direct Urea process (2), there is a great nec¬ 
essariness for an accurate measurement of the neutron skin. 
The information about the neutron skin thickness of 208 Pb 
A R np = 0.156^qq2i fm from proton inelastic scattering (33) , 
A R np = 0.17 ± 0.03 fm from chiral effective field theory 11341 
and A R np = 0.191 + 0.032 fm from a -decay energies (35) 
shows a big challenge of a precise constraint about the neu¬ 
tron skin thickness. 


2 Method 

For heavy nuclei, such as 208 Pb, the nuclear surface region 
where nucleon density is much less than the saturation den¬ 
sity, contributes dominantly to symmetry energy {35) ■ So it 
is better to describe the density-dependent symmetry energy 
by expanding around the density that is below the satura¬ 
tion density. Around the nuclear matter subsaturation den¬ 
sity p = 0.1 lfm \ the symmetry energy S ip) is expanded to 
second order in term of density p as 


Sip) - 5o.ii 


Tq.ii / P - 0.11 \ 
+ 3 \ 0.11 ) 


K sym ip- o.n \ 2 

+ 18 \ 0.11 ) ' 


( 1 ) 


There is an available connection that the a sym iA) of finite nu¬ 
clei is approximately equal to S ip a) of the nuclear matter at 
a reference density pa, which is proposed by Centelles, et 
al (30). It makes the symmetry energy of the nuclear matter 
in contact with the one of finite nuclei, and thus allows one 
to explore the density dependence of the symmetry energy 
S ip). We have extracted the mass-dependent symmetry en¬ 
ergy coefficients a sym iA) in our previous work and found that 
asymiA) — 22.4 MeV for 208 Pb (36) . The reference density 
p A for 208 Pb is figured out to be 0.55po = 0.088 fm -3 with 
po = 0.16 fm~ 3 . Substituting the pa in Eq.[l], we obtain 


22.4 MeV = S 0 .„ 


0.6T 0 I i - 0.02K iym 

9 


( 2 ) 


There are several shapes of symmetry energy as a function of 
density (37) . It is found that the expression Sip) — S o(p/po) r 
or Sip) = 12.5(p/po) 2/3 + Cpip/po) y is not universal. But 


the shape from the density-dependent M3Y interaction {38) 
is much better. With it, we can describe the symmetry energy 
S ip) in nuclear matter as the following formulism 

S(p> = Ck (oir)’ + Cl (out) + Cl (oil)’ • (3) 

where Ck = h ^' n - 9.6 MeV, Ci and C 2 are two param¬ 

eters that are required to be determined. Lo.n = 3p|^|o.n 
and K sym — 9p 2 |p-|o.n are slope and curvature parameters 
at the density of 0.11 fm 2 respectively. There is an use¬ 
ful connection among the symmetry energy at the density of 
0.1 lfnT 3 So.n, Lo.n and K sym 


K sy m - 5Lq.11 - 15So.ii + 28.8, (4) 


which is sufficient to estimate the small contribution of K sym 
in Eq.[l]. As a consequence, Eq.[l] can be rewritten as 

, 195.24 -0.52L 0 .n L 0 . u Ip - 0.11 

5(p)= —9— + —hur- 

4.14Lo. u -317.79 /p - 0.11 \ 2 
+ 18 \ 0.11 ) ' 

Therefore, once the Lo.n is determined, the symmetry energy 
Sip) at subsaturation density is obtained. Thus, the centre 
goal of this work is to determine Lo n- In the local density 
approximation, the symmetry energy coefficient a sym iA) can 
be calculated (39) as 

a sy m(A) = -^ f d\pir)S[pir))[5ir)}\ (6) 

where Xq — (V - Z)/A is the whole nuclear isospin asymme¬ 
try, <5(r) = ip n ir) - Ppir))/pir) the isospin asymmetry profile, 
pir) the whole nucleon density and a sym = 22.4 MeV for 
208 Pb. p„(r) and p p ir) are the density profile of neutron and 
proton respectively in nucleus that can be given by the mean 
field models such as Skyrme-Hartree-Fock. Here, we calcu¬ 
lated pir), Sir) and the neutron skin thickness A R np with the 
Skyrme interactions |40l : Sill, SLy4, SLy5, SkI2, SkI4, SkP, 
SkM*, SGII, T4, T6, BSk8, LNS1, LNS5, HFB17, KDE0, 
KDE, SKz2, SKz4, SV, MSkl, MSkA, v090, SK255, SK272, 
SIV, SkMP, MSL0, SKA, SKSC15 and SKSC40. Substitut¬ 
ing Eq.[5] in Eq.[6], the slope parameter Lo.n is determined 
and the obtained Lo.n with this method is labeled L^’j { . In¬ 
deed, in order to extract the nuclear symmetry energy that re¬ 
lates solely to the nuclear force, the effect due to the fact that 
the Coulomb interaction effectively polarizes the neutron and 
proton densities should be subtracted. So we calculate the 
neutron density p,,ir) and proton density p p ir) by switching 
off the Coulomb interaction. On the contrary, the Coulomb 
interaction is necessary to work out the neutron skin thick¬ 
ness A R np on account of that the protons and neutrons move 
slightly away from the core as a result of the Coulomb repul¬ 
sion between charged protons. 
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3 Results 

Known that a sym (A) — 22.4 of heavy spherical nucleus 208 Pb 
above, we can obtain Lo n with various Skyrme interactions. 
The results given by some typical Skyrme interactions men¬ 
tioned above are shown in Table 1. Ly 6 ^ and j are the 
slope parameter and the symmetry energy at the density of 
0.1 lfnT 3 calculated with the properties of finite nuclei while 
Lo.ii and So.n are those calculated with the Skyrme interac¬ 
tions for nuclear matter directly. Lo.i i obtained by the Skyrme 
procedure for nuclear matter directly has a large scope but 
the L® j computed with Eq.[6] ranges narrowly. And it is 
the same for the comparison between both symmetry energy, 
namely S ( ( , 6 [ j and .S’o.i i gained in the two ways. It thus appears 
that the method of calculating the features of nuclear matter 
by using the structure properties of the finite nuclei is feasible 
and effective. 

Table 1 Comparison between the values of the slope parameter and the 
symmetry energy at the density of 0.1 lfm -3 with the properties of finite nu¬ 
clei namely, Eq.[6], and those calculated with Skyrme interactions directly. 


Force 

r ( 6) 

^0.11 

o (6) 

13 0.11 

To. n 

5 0.11 


(MeV) 

(MeV) 

(MeV) 

(MeV) 

sin 

50.71 

25.54 

32.02 

26.06 

SLy4 

48.03 

25.34 

42.09 

26.49 

SLy5 

47.29 

25.28 

43.52 

26.22 

SkI2 

60.30 

26.27 

68.17 

23.23 

SkI4 

49.74 

25.47 

46.09 

22.9 

SkP 

46.33 

25.21 

35.26 

26.21 

SkM* 

50.33 

25.51 

44.17 

24.32 

SkMP 

53.31 

25.74 

53.45 

22.59 

soil 

45.26 

25.13 

37.95 

22.18 

T4 

59.63 

26.22 

67.33 

25.63 

T6 

52.77 

26.70 

38.58 

25.43 

BSk8 

50.85 

25.55 

28.82 

25.18 

LNS1 

44.54 

25.07 

38.54 

25.31 

LNS5 

45.17 

25.12 

44.91 

23.14 

HFB17 

49.00 

25.41 

40.08 

25.25 

KDE0 

46.42 

25.22 

43.89 

27.30 

KDE 

44.87 

25.10 

41.05 

26.39 

SKz2 

43.53 

25.00 

33.36 

28.69 

SKz4 

43.15 

24.97 

25.59 

29.90 

SV 

42.81 

24.94 

66.72 

23.62 

MSkl 

49.15 

25.42 

39.69 

25.49 

MSkA 

50.07 

25.49 

50.90 

24.31 

v090 

43.35 

25.06 

27.44 

25.84 

SK255 

58.60 

26.14 

71.00 

27.56 

SK272 

58.72 

26.15 

70.15 

28.13 

SIV 

52.04 

25.64 

55.59 

24.91 

MSL0 

52.29 

25.66 

49.38 

23.15 

SKA 

55.55 

25.91 

58.98 

25.25 

SKSC15 

43.15 

24.97 

27.99 

25.58 

SKSC40 

40.51 

24.77 

21.12 

25.98 


All the results about Lo.n from Eq.[6] (labeled L®) are 
shown in the left panel of Fig. 1 between the two dashed lines. 
The calculated Lo.n directly for the nuclear matter is shown 
in the horizontal axis and the the neutron skin thickness A R np 
is shown in the vertical axis. The linear relation is demon¬ 
strated obviously. By employing the least square fitting, we 
obtained the correlation between AR„ P in 208 Pb and Lo.n as 
following 

A R„ p = (0.0279 ± 0.00338) + (0.00307 + 0.00007)L 0 . n (7) 

where A R np and Lo.n are in units of fm and MeV respec¬ 
tively. The gained linear fitting coloured red is fine with the 
correlation coefficient r up to 0.984. The slope paramter of 
symmetry energy at the density of 0.1 lfnT 3 L[j n given by 
our present work with the finite nuclear properties ranges 
from 40.5 MeV to 60.3 MeV within the 30 sets of Skyrme 
interaction. Consequently, the neutron skin thickness A R np 
of 208 Pb can be predicted in terms of this correlation, which 
is constrained from 0.15 fm to 0.21 fm. If A R np of 208 Pb is 
greater than 0.24 fm, the direct URCA process to cool down a 
1.4M 0 neutron star is allowed ESI . Our result is too small to 
make the direct Urea process in the 1.4M 0 neutron star occur. 
There are some recent data about the neutron skin thickness 
AR np of 208 Pb, which are shown in Table 2, obtained from 
various approaches. It appears that our result consists well 
with 0.180 + 0.035 fm from the pygmy dipole resonances 
(PDR) m and 0.185 + 0.03 fm from mean field Q7). Re¬ 
cently, Dong et al employed a more practicable strategy than 
the current lead radius experiment (PREX) to probe the neu¬ 
tron skin thickness of 208 Pb based on a high linear correlation 
between the AR np and .1 - a synl m, where J is the symmetry 
energy (coefficient) of nuclear matter at saturation density. 
The obtained A R np in 208 Pb was 0.176 + 0.021 fm robustly, 
being consistent with the present result. 

Table 2 The neutron skin thickness AR„ p of 208 Pb probed in various inde- 
pendent studies. 


Reference Method AR np (fm) 


(SI 

proton inelastic scattering 

0 156 +0025 
u. iju_o 021 

ED 

proton elastic scattering 

0 211 +0 054 

U -0.063 

E) 

chiral effective field theory 

0.17 ±0.03 

EU 

alpha-decay energies 

0.191 ±0.032 

(SI 

mean field 

0.185 ±0.035 

ED 

pygmy dipole resonance 

0.185 ±0.035 

ED 

pygmy dipole resonance 

0.194 ±0.024 

Present 

Nuclear mass differences 

0.15-0.21 


In addition, we also show the correlations of the neutron 
skin thickness A R„ p of 208 Pb and the slope coefficient L () of 
symmetry energy at the saturation density in the right panel 
of Fig. 1 with a linear fitting. The linear relation is strong with 
the correlation coefficient r=0.964 but is weaker than the one 
in the left panel. It can be concluded that the neutron skin 
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thickness of heavy nuclei is uniquely fixed by the symmetry 
energy density slope at a subsaturation density of 0.1 lftrT 3 
rather than at the saturation density, which is in agreement 
with the conclusion in |45| . 

After we have got the the slope coefficients Ll 6 j ; of symme¬ 
try energy at the density of 0.11 fm -3 by 30 sets of Skyrme pa¬ 
rameters, with Eq.[2] and Eq.[4], the symmetry energy 


and the curvature paramter K sym at the density of 0.11 fm 3 
can be measured. The calculated is 24.8-26.3 MeV 
and /f sym is from -150.1 MeV to -68.3 MeV. Then, the de¬ 
scription Eq.[l] for the density-dependent symmetry energy 
around the subsaturation density of 0.1 lfm~ 3 becomes clear. 
The behaviors of established symmetry energy S (p) are pre¬ 
sented as a function of density p in Fig. 2. 




Figure 1 Correlation of the neutron skin thickness A R np of 208 Pb and the slope parameter Lo n (Lq) of symmetry energy at the density of 0.1 lfm 3 (saturation 


density) in the left (right) panel within the Skyrme interactions. 

4 Summary 


In summary, we have employed the symmetry energy coeffi¬ 
cient a sym (A) extracted by experimental nuclear masses in our 
previous work and the density profiles of heavy nuclei to ex¬ 
plore the behavior of the density-dependent symmetry energy 
around the subsaturation density of 0.1 lfnV 3 . The estimated 
value of the slope parameter Lq i i of symmetry energy at the 
density of 0.1 lfnT 3 is 40.5-60.3 MeV. With the nice linear 
correlation of Lo n and the neutron skin A R np of 208 Pb , A R np 
is predicted to be 0.15-0.21 fm, which is too small to allow 
the direct Urea process in the 1.4M G neutron stars. It proves 
that the approach of probing the nuclear matter features by 
using the information about the nuclear structure is effective. 
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